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Abstract 

The  bounded  interval  generalized  assignment  problem  is  exemplified 
by  the  problem  of  assigning  tasks  to  agents  so  that  each  task  is  assigned 
to  exactly  one  agent  and  the  time  required  to  complete  the  set  of  tasks 
assigned  to  any  one  agent  falls  between  prespecified  lower  and  upper  bounds. 
This  note  describes  an  efficient  algorithm  for  solving  this  problem. 


1.0  Introduction 


Assignment  models  are  characterized  by  a set  of  indivisible  tasks,  a set 


of  agents  responsible  for  completing  the  tasks,  and  a criterion  function  for 


evaluating  pairings  of  agents  and  tasks.  In  many  practical  problems,  the 


number  of  tasks  exceeds  the  number  of  agents,  and  this  characteristic  gives 


rise  to  "many-for-one"  assignment  models.  In  these  models,  it  is  possible  for 


an  agent  to  be  assigned  a set  of  tasks  which  collectively  require  more  of  some 


resource  (e.g.  time)  than  the  agent  has  available.  To  accommodate  this  factor 


assignment  models  and  special  purpose  algorithms  have  been  proposed  that 


explicitly  consider  the  resource  burden  associated  with  multiple  task  assign 


ments  [ 4,ll].  An  important  generalization  of  these  models,  known  as  the 


generalized  assignment  problem,  additionally  considers  differences  in  the  agents 


abilities  to  complete  the  same  task  [ 9].  Specific  applications  of  this  model 


include  the  allocation  of  inventory  items  to  warehouses  [ 4 ] , the  assignment  of 


ships  to  shipyards  for  overhaul  [8],  the  assignment  of  jobs  to  computers  in 


computer  networks  [ 1]  and  the  assignment  of  software  development  tasks  to 


Analysis  of  other  assignment  problems  reveals  additional  restrictions  limit 


ing  the  pairings  of  agents  and  tasks.  As  suggested  in  [12],  some  situations 


require  each  agent  to  contribute  a minimal  amount  of  his  resources  to  completing 


the  tasks.  Such  constraints  may  reflect  managerial  policies  that  require  an 


equitable  distribution  of  work.  Consideration  must  be  given  to  the  relative 


strengths  and  weaknesses  of  the  agents,  and  the  agents  should  neither  be  over 


burdened  nor  underutilized  [ 5] 


Analagous  restrictions  arise  in  other  contexts.  In  machine  loading  models. 


be  desirable  to  balance  the  machine  workloads  rather  than  allowing  some 


heavily  loaded  and  some  lightly  loaded  machines.  Similarly,  capacity  and 
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configuration  constraints  for  facility  location  models  may  restrict  both  the 
minimum  and  the  maximum  size  of  a facility  or  the  number  of  facilities  [10]. 
Finally,  territory  design  procedures  for  problems  of  political  districting, 
school  districting  and  sales  districting  [14]  require  an  equitable  distribution 
of  some  entity  (such  as  voters,  minority  students  or  sales  potential)  among 


the  districts. 

The  purpose  of  this  note  is  to  formulate  these  generalized  assignment 
problems  with  bounded  interval  resource  limitations  and  to  show  how  the  branch 
and  bound  algorithm  of  [ 9]  may  be  extended  to  solve  these  problems. 


2 . 0 The  Bounded  Interval  Generalized  Assignment  Problem 

The  bounded  interval  generalized  assignment  problem  may  be  formulated  as: 


(P)  minimize 


ZEE  Ci1Xij 
iel  jeJ  3 3 


subject  to  E x..  = 1 for  every  jeJ 
iel  1J 


a.  £ E r x.,  <_  b for  every  iel  (2) 

1 jeJ  1J  1J 

x^j  = 0 or  1 for  every  iel,  jeJ  (3) 

The  set  J:  { 1 , 2,  . . . , n}  represents  the  tasks  which  must  be  completed,  the 
set  I = {1,  2,  ...,  m}  represents  the  agents  who  can  be  assigned,  and  the  parameters 
c„  and  r^j  are,  respectively,  the  cost  incurred  and  the  resource  required  if 
agent  i is  assigned  task  j.  The  natural  interpretation  of  the  decision  variable 


r 1 of  agent  i is  assigned  task  j 
xij  0 otherwise 

The  constraints  (1)  and  (3)  insure  that  each  task  is  uniquely  assigned  to  a 
single  agent.  Constraints  (2)  insure  that  each  agent  expends  no  less  than  a^ 
and  no  more  than  b^  resource  units  in  accomplishing  his  assigned  tasks.  The 
usual  generalized  assignment  problem  [ 9]  is  a special  case  of  (P)  in  which 
= 0 for  all  iel.  We  shall  refer  to  that  problem  as  (PQ) . 
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3.0  An  Algorithm  For  Solving  ( P ) 

The  branch  and  bound  algorithm  of  [ 9]  exploits  the  special  mathematical 
structure  of  (Pq).  Our  objective  here  is  to  show  how  that  algorithm  can  be 
easily  extended  to  solve  (P) . Moreoever,  we  shall  show  how  the  necessary  modi- 
fications preserve,  as  much  as  possible,  the  computational  efficiency  of  the 
Ross  and  Soland  algorithm. 

A bounding  procedure  similar  to  that  of  [ 9]  is  employed.  Its  validity 
can  be  established  using  the  theory  of  Lagrangean  relaxation  [ 6],  That  is, 
lower  bounds  for  (P)  can  be  calculated  by  solving  the  following  relaxation  of 
(P): 

(PR, ) minimize  z.  = Z Z c x + Z A (1  - Z x. .) 

A A iel  jeJ  i3  ij  jeJ  j iel 


subject  to: 


a < Z r , x _<  b for  every  iel 

±_  jeJ  3 3 


Xij 


0 or  1 for  every  iel,  jeJ. 

The  problem  (PR^)  relaxes  the  requirement  that  every  task  must  be  completed,  but 
it  incorporates  the  constraints  (1)  in  the  objective  function  with  fixed  multi- 
pliers to  penalize  the  optimal  value  whenever  these  constraints  are  not  satisfied. 
The  objective  function  of  (PR^)  may  also  be  written  as 

Z\  ~ E A - maximum  [ Z Z x ]. 

A jeJ  3 iel  jeJ  3 13  13 

With  the  objective  function  in  this  form,  it  is  clear  that  (PR^)  separates  into  a 

collection  of  "bounded  interval"  binary  knapsack  problems,  one  for  each  iel.  Each 

of  these  knapsack  problems  has  one  constraint  of  the  form  (2) . Such  constraints 

differ  from  the  usual  knapsack  constraint  in  that  the  summation  is  bounded  both 

from  below  and  from  above.  Thus,  the  efficiency  of  the  bounding  procedure  rests 

in  part  on  the  fact  that  bounded  interval  binary  knapsack  problems  can  be  solved 

very  efficiently  ( 3,  13]. 
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As  in  [ 9],  each  X in  the  objective  function  of  (PR.)  is  set  equal  to  c_ . , 

J A Z j 

the  second  smallest  value  of  c^  for  all  iel.  These  X^  are,  in  turn,  optimal 

dual  multipliers  for  the  trivial  bounded  variable  linear  program: 

(PR^)  minimize  Z^  = E E c^.x 

iel  jeJ  ^ ^ 

subject  to:  E x..  = 1 for  every  jeJ 

jeJ  1J 

0 <_  x <_  1 for  every  iel,  jeJ. 

Other  values  for  the  X^  could  be  used,  but  the  ones  just  described  have  the  ad- 
vantage of  being  easy  to  compute. 

Thus,  determining  a lower  bound  requires  two  steps.  First,  solve  (PR^) , then 

solve  (PR^).  If,  however,  the  primal  solution  X = (x^ ) to  (PR^) , which  satisfies 

(1),  should  also  satisfy  (2),  then  Z = = Z^,  and  (PRp  need  not  be  solved. 

Ordinarily,  X will  not  satisfy  (2),  and  (PR^)  must  be  solved  to  obtain  Z^. 

In  addition  to  providing  a lower  bound,  the  solutions  to  (PR^)  and  (PR^)  are 

used  in  selecting  a branching  (or  separation)  variable  for  defining  subsequent 

candidate  problems  (CP).  As  noted  above,  the  solution  to  (PR^),  X,  is  usually 

not  feasible  to  (2).  In  essence,  the  solution  to  (PR.),  X = (x  ) , may  be  inter- 

x lj 

preted  as  recommending  changes  in  X which  must  be  made  in  order  to  satisfy  (2) . 

That  is,  it  is  possible  that  for  some  jeJ,  E x = 0 to  avoid  overloading 
_ iEl  1J 

any  agent  or  E x,  > 1 to  insure  every  agent  uses  a minimum  amount  of  his  re- 
iel  = 

source.  Those  variables  x^  with  an  optimal  value  of  one  indicate  agent-task 
pairings  that  should  be  made  whereas  those  x^  with  an  optimal  value  of  zero 
indicate  pairings  that  should  be  avoided.  Thus,  these  variable  values  indicate 
changes  that  will  reduce  the  aggregate  infeasibility  of  x in  (2),  and  they  are 
helpful  in  choosing  a branching  variable. 
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To  formalize  the  concept  of  reducing  aggregate  infeasibility,  we  define  the 
infeasibility  in  constraint  i prior  to  taking  a branch  to  be 
= max  {0,  d+,  d^ 

where  d+  - ^ -b±, 

di  = ai  “Jj  rijXi.1’ 

The  set  I+  ={iel|d*  >0}  identifies  those  constraints  (2)  for  which  X exceeds  the 
upper  bound,  and  I = (iel|d^  >0}  identifies  those  constraints  (2)  for  which  X 
fails  to  satisfy  the  lower  bound. 

Suppose  I+  t 0 and  kefjejlx^.  = 1 and  iel+};  if  x^  Is  set  t0  0 then  d+  and 


1 


d^  become: 


di  1 rijxij  " bi  " rik 

J CJ 


di  ■ ai  - ru  xu  + rik 


and  the  resulting  infeasibility  in  constraint  i becomes 


= max  {0,d+,  d^. 


Assuming  that  task  k is  reassigned  to  the  second  least  costly  agent,  (say  agent  h, 

min 

1*1 


where  c = min  c,,  ) then  the  infeasibility  in  constraint  h becomes 
hk  ...  Ik 


= max  {0,d^,dh> 


where 


du  = Z rv.  4 XM 

h jeJ  hj  hj 


K + r.  , 

h hk 


dh  " ah  ~JeJ  rhj  Xhj  " rhk. 


Hence  the  net  difference  in  total  infeasibility  is: 

4DJ " (Di  + V - <Di  + DS> 

If  AD^  > 0 then  setting  x^  = 0 yields  a reduction  in  aggregate  infeasibility, 
and  if  AD^  < 0 then  such  a branch  will  not  reduce  aggregate  infeasibility. 

Similarly,  suppose  that  I ^0  and  ke{jej  | x = 0 and  iel  };  if  x^  were 
set  to  1 then  d^  and  d^  become: 

di  ru  ' bi + rik 

di  ' “i  'jjj1!!  *ij  “ rlk 

and  the  resulting  infeasibility  in  constraint  i would  be 

= max  {0,  d*,  d^}. 


If  x^  is  set  to  1 then  task  k is  assigned  to  agent  i and  agent  i^  relinquishes 
it,  where  i^  = min  c^-  Hence  the  infeasibility  for  constraint  i^  becomes 

k 4*  — 

D = max  {0,  dj_  , d.} 

k k xk 


where 


dt  = E rj  iXj  b.  - r 

k j£j  V V \ ‘kk 


V \ ■ iL  rw+  v • 


The  net  difference  in  infeasibility  is 

ADi  " <Di  + Di  > " (Di  + > 

k Xk 

where  D±  and  Dlfc  are  the  infeasibilities  in  constraints  i and  ifc  prior  to  any 

i, 

branch.  As  before,  if  ADj^  > 0 then  there  is  reduced  infeasibility  following 
a branch  on  variable  x,,  . 


I 
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Several  rules  for  selecting  the  branching  variable,  x.*,*,  are  formulated 


as  follows: 


I.  a)  x^*j*  is  that  variable  for  which 

AD^  = max  {AD^} 

1 (i,  j)  e H+UH'  1 

where  H+  = { (i , J ) I x£ j = ® ant^  iEl+} 

H"  = {(i,  J)|x  = 1 and  iel-} 

b)  If  AD^  = 0 in  a)  then  xi*j*  is  that  variable  for  which 


where 


(i,  j)  £ (G  -H  ) U (G"-H_) 
G+  = {(i,  J) |*  = 1 and  iel+} 


{AD^} 


G"  = ((i,  j) 


0 and  iel  }. 


II.  a)  x. ...  is  that  variable  for  which 
1 J r 


Pi*,*  = mln 


(i,  j)eG+  [ Ad| 


Ycu 


(i,  j)eG  |_  AD£ 


Cij"Cikj 


b)  If  Ad|  = 0 for  all  (i,  j)eG+UG  then  x^*j*  is  that  variable  for  which 


(i,  j ) £E 


* V 


where  E = {(i,  j)|x  = 1 and  iel  }, 

denotes  the  set  of  tasks  assigned  to  agent  i by  prior 
branching. 

Rules  la  and  lb  are  designed  to  choose  that  variable  which  reduces  the  post  branch 
aggregate  infeasibility  by  the  greatest  amount.  Rule  Ila  conditions  the  choice  of 
branching  variable  on  the  additional  cost  incurred  per  unit  reduction  in  infeasi- 
bility. Rule  lib  is  the  one  used  in  [ 9];  the  variable  chosen  by  this  rule  repre- 
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sents  an  agent-task  pairing  which  should  be  made  considering  the  penalty  for  not 
doing  so  weighted  by  the  fraction  of  the  agent's  remaining  free  resources  con- 
sumed by  the  assignment. 

As  the  algorithm  progresses  and  new  candidate  problems  (CP)  are  defined 
by  the  branching  process,  the  additional  steps  given  below  may  be  taken  to 
facilitate  fathoming.  These  steps  are  specialized  adaptations  of  more  general 
forcing  (or  variable  fixing)  tests  suggested  by  Balas  [ 2]  and  Glover  [ 7]. 


In  solving  any  (CP),  any  x^,.,  for  which  r^,.,  > b^,  - E xi'iri'-  may 

3 1 J 1 jeFit  1 J 1 2 

be  set  equal  to  zero.  Here  F^,  denotes  those  jeJ  for  which  x^,j  has  been  assigned 

a value  of  zero  or  one  by  prior  branching  or  variable  fixing  tests.  Similarly, 

if  there  is  an  x . , , , for  which  a.,  - E x, ,.r. ,.  > E r. ,.-r. 

1 J 1 1 * 1 i JEJ-Fj,  1 J J 

then  x^,^,  must  be  set  equal  to  one  in  the  solution  to  (CP).  These  variable 
forcing  tests  may  subsequently  result  in  other  variables  being  forced  to  zero 


or  to  one  when  all  of  the  resultant  implications  are  considered.  Moreover, 


forcing  certain  variables  to  zero  or  to  one  in  the  solution  to  (CP)  may  affect 

I 

the  values  of  some  of  the  A^  obtained  from  solving  (CPR^) . This  change  may, 

in  turn,  increase  the  value  of  the  lower  bound  provided  by  (CPR^) . 

Another  test  may  be  used  to  check  the  feasibility  of  (P)  (or  any  candidate 

subproblem).  Summing  the  constraints  (2)  together  yields  the  constraint  (5): 

A = E a.  £ E E r x < E b,  = B (5) 
iel  iel  jeJ  J J iel 

This  new  constraint,  together  with  constraints  (1),  implies  that  for  any  feasible 

soltuion  to  (P) , we  must  have: 

E r'  _>  A and  E r",  B (6) 

jeJ  J jeJ  J 


where 
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The  values  necessary  for  the  tests  (6)  can  be  updated  easily  as  part  of  the 
branching  process  in  order  to  apply  this  test  to  each  (CP). 

The  algorithm  terminates  in  the  usual  way  when  all  candidate  problems  have 
been  fathomed. 


4 . 0 Conclusion 

This  note  has  described  an  efficient  branch  and  bound  algorithm  for  the 
bounded  interval  generalized  assignment  problem.  The  algorithm  serves  as  a 
useful  tool  for  solving  a large  number  of  applications  of  this  assignment 
model,  a representative  sample  of  which  are  mentioned  in  the  introduction. 
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